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ABSTRACT 



The observed distributions of the thermodynamic variables near stratocumulus 
top are highly bimodal. We examine two simple models of sub-grid fractional 
cloudiness motivated by this observed blmodallty. In both models, certain low- 
order moments of two independent, moist-conservative thermodynamic variables are 
assumed to be known. 

The first model is based on the assumption of two discrete populations of 
parcels: a warm-dry population and a cool-moist population. If only the first 

and second moments are assumed to be known, the number of unknowns exceeds the 
number of independent equations. If the third moments are assumed to be known 
as well, the number of independent equations exceeds the number of unknowns. 

The second model is based on the assumption of a continuous joint bimodal 
distribution of parcels, obtained as the weighted sura of two binorraal distribu- 
tions. For this model, the third moments are used to obtain 9 Independent non- 
linear algebraic equations in 11 unknowns. Two additional equations are needed 
to determine the covariances within the two subpopulations. In case these two 
internal covariances vanish, the system of equations can be solved analytically. 



1 . Introduction 


Stratocumulus sheets are one of the most widespread cloud forms over the 
globe. For this reason, and In view of the Importance of cloudiness as a fore- 
cast product for numerical weather prediction, and as a pivotal component of 
the climate system (JOC, 1975), major efforts have been mounted to simulate 
stratocumulus clouds In general circulation models (Randall, 1976, 1982 j Sllngo, 
1980; Ramanathan and Dickinson, 1981; Suarez et» al. , 1983). These efforts have 
been limited by a number of simplifying assumptions, notably Including the 
assumption that at a given level the fractional cloudiness Is either zero or 
one. 

Fractional cloudiness enters the problem In two different ways. First, as 
discussed by Sundqvlst (1978, 1981), cloud decks are rarely unbroken over an 
area the size of a general circulation model grid box (several hundred kilometers 
on a side). Sundqvlst has proposed methods to determine the macroscale fractional 
cloudiness within a grid box, and this will undoubtedly be an active area of 
research In coming years. 

However, fractional cloudiness also enters Into the micro-scale dynamics 
governing the local structure of the clouds. It Is this aspect of the frac- 
tional cloudiness problem that Is addressed in the current proposal. Eradiation 
and the buoyancy flux, which strongly control the rate of cloud-top entrainment 
(Randall, 1980a, b) , depend sensitively on the fractional cloudiness at each 
level. Near cloud top, there Is a layer within which the fractional cloudiness 
decreases with height from near 100% to zero. 

The upper boundary of a layer of stratocumulus clouds Is Inevitably somewhat 
Irregular; turrets can be pushed up here and there by convective turbulence, and 
shearing Instabilities can give rise to wave-llke features. The bumpiness of 


cloud top complicates the theoretical description of such Important cloud-top 
processes as radiative cooling and turbulent entrainment, since at a given level 
the fractional cloudiness is greater than zero but less than one* The clear air 
near cloud top is typically very warm and dry compared to the cool, moist cloudy 
air lying nearby at the same height (Fig. 1). 

Up to now, only the high-resolution three-dimensional model of Deardorff 
(1980) has explicitly simulated the effects of partial cloudiness near cloud 
top; Deardorff remarked that even in his model the vertical resolution (50 m) 
appeared to be Inadequate to accurately represent processes near cloud-top* He 
advocated the use of a parameterization of sub-grid fractional cloudiness as an 
aid in the simulation of cloud-top processes* 

Such parameterizatlons have been proposed and discussed by Deardorff and 
Someria (1977), Mellor (1977), Oliver et* al* (1978), and Bougeault (1981a, b)* • 
The first two studies dealt with fractional cloudiness in air characterized by 
binormal probability distributions of a pair of moist-conservative thermodynamic 
variables, such as liquid-water static energy and total mixing ratio. Oliver 
et* al* proposed a simpler model, based on ad hoc assumptions* Bougeault 
generalized the Deardorff-Someria-Mellor model to allow skewed distributions, 
at the expense of additional tunable parameters* 

However, the observed distributions of the thermod 3 mamic variables near 
stratocumulus cloud-top are not well-represented by gausslan or even skewed- 
gaussian distributions, because of the Very sharp changes of temperature and 
moisture with height in the capping inversion (Lumley, 1979). The observed 
distributions are stongly blmodal; there are a lot of warm-dry parcels from the 
inversion above, and a lot of cool-moist parcels from the cloud below, but 
relatively few warm-moist or cool-dry parcels. An example given by Mahrt and 


Paumiec (1982) Is shown in Fig. 2. A parameterization of the subgrid fractional 
cloudiness near stratocumulus top must take this blmodality into account. 

The purpose of this paper is to report on efforts to develop a new para- 
meterization of subgrid-fractional cloudiness, in which joint blmodality of 
the thermodynamic variables Is permitted. The parameterization Is Intended 
to represent the subgrld-f ractlonal cloudiness near the top of a stratocumulus 
deck; however, It may have wider applications. The results of this study should 
make ‘possible Improved parameterlzatlons of stratiform cloudiness for use In 
climate models. They may also provide some guidance for parameterlzatlons of 
macro-scale fractional cloudiness for climate models. 



OHlGlNAt 
OF POOR 

We begin by investigating a highly Idealized model for which blmodallty 
Is the whole story. Consider a mass of air consisting entirely of two kinds 
of parcels. Parcels of Type 1 have liquid water static energy sj and total 
mixing ratio (vapor plus liquid) qj; parcels of Type 2 have liquid water 
static energy 33 nnd total mixing ratio q 2 * Since (s,q) can assume only 
the discrete values (sj^, q^) and (s 2 , <\ 2 ) * refer to this as the "discrete 
model." At a given height, the fraction of Type 1 air Is f^, and the fraction 
of Type 2 air is (1 - f^). The physical picture is illustrated in Fig. 3. 

The joint probability density function Is zero except at a pair of points, 
where It Is infinite (Fig. 4). 

The first and second moments of s and q are given by 

s = s^f^ + s^(l - f^) , ( 2 . 1 ) 

q = q^f^ + ~ ’ (2.2) 

s'T = (s^ - 7)2,f^ + (s^ - 1)2(1 - fj) . (2.3) 

q’2 = (q^ - q)2f^ + (q^ - q)2(l - f^) , (2.4) 

s'q' * (s^ - s)(q^ - q)fj^ + (s^ s)(q 2 “ q)(l “ f^). (2.5) 

We assume that these moments are known, as for example In a seconu-order 
closure model. Then (2. 1-2,5) appear to represent 5 equations which c«n be 
solved for the 5 unknowns s^, S 2 , qi, q 2 f fl. 
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From (2. 1-2.4) we can solve for * 2 * *l2 terms of f^, as 


follows: From (2.1), we obtain 


- s - (s^ - Sj,)(l - f^) , 


(2.6) 

-s.(s^-Sj)£, . 

ORiGINAL PAGE IS 
OF POOR QUALITY 

(2.7) 


Substituting into (2.3), we find that 

(s^ - s^)?' - - fj)]. (2.8) 

Going back to (2.6), we square both sides and substitute from (2.8) to 
obtain 

8^ - ? ± {[(1 - f^)/f^l s~’2)l''2 . (2.9) 

Similarly, from (2.7) and (2.8), we find that 

s^ - ? ± Itf^/d - f^)l . (2.10) 

By analogy, and q£ are given by' 

q^ - q ± ([1 - f^)/f^I ^}l/2 ^ (2.U) 

92 " ^ * {[fj/(l ” fj)] q’2}l/2 . (2.12) 


If we choose the -f in (2.9), then sj > s necessarily, and so S 2 must be less 
than T, i.e. , we must choose the minus In (2.10). A similar argument holds for 
(2.11) and (2.12). Suppose we choose the (+) in (2.9), so that sj > s. Then 
from (2.5) we see that choice of the (+) in (2.4), implying qj > ^» would 
guarantee s'q' > 0. Thus if the sign of s'q* is known, the choice of signs in 
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(2.9) and (2.11) cannot be made independently. Then let's proceed as follows. 
Define 


1, S'q*' > 0 

OF POOR 

M = 

- 1 , < 0 . 

We arbitrarily choose the + in (2.9), and all else follows: 

= 7 + {[(1 - s ' 2 } i /2 ^ 

=7- - f^)] s'2}l/2 , 

«= q + M {[(1 - f^)/f J q'2}l/2 ^ 

q^ = q - M {[f^/(l “ f^)] q'2}i/2 


(2.13) 

(2.14) 

(2.15) 

(2.16) 
(2.17) 


Note that now by definition Type 1 air has a higher liquid water static energy 
than Type 2 air. 

From (2.14-2.17) we can obtain sj, S 2 » qx, and q 2 if fi is known. Thus 
we rely on (2.5) to determine fj. Substituting (2.14-2.17) into (2.5), we 
obtain 

7'V' = M»(s'2q'^)^-^2 , (2.18) 

Since fj has dropped out of (2.18), knowledge of the covariance does not help 
us to determine fj. We have more unknowns than we have equations. Some new 
information is needed. 

One possibility is to assume that- the predicted variances are explained 
by the smallest possible differences between the two types of air, i.e. assume 
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that <8^ - minimized. However, this leads to ■ 1/2 in 

every case, so the' assumption is not useful. 

Another possibility is, to Introduce the third moments: q'^, s'^q', 

and s'q'^. These we assume to be known, as for example in a third**order 
closure model. We find that 


s»3 , - s)3f^ + (s^ - s)^ (1 - fj^) 

- {((1 - f, - {l£/(l - £j)l - fj). 


(2.19) 

From (2.19) we can derive a quadratic equation for fp 

fj2[4 + (7^^)2(7T2)-3] . + (I^)2(7T2)-3] + 1 - Q; (2.20) 

the solution is 

t 

f ml {[1 -7^ [4(T^)3 + (7^)2 1-1/2} , (2.21) 

2 

In (2.21), the sign before the radical has been chosen by the following consi- 
derations. Regardless of the choice of sign, we have fj - 1/2 for s' ^ ■ 0. 

Fig. 5 shows a section along which fj < 1/2. Suppose for simplicity that sj^ > 0 
and S 2 * “Si. (This is consistent with our choice sj > S 2 .) It is clear from 
the figure that s'^ > 0, since along the section there la a small amount of air 
with large positive s', and a large amount with small negative s'. We conclude 
that fj < 1/2 for s'^ > 0, and fj^ < 1/2 for s'^ > 0. This explains the 
choice of the minus sign in (2.21). 

Since we need 0 < f^ ^ 1, we must require that 


-1 < s'3 [4 (s' 2)3 + (gi3)2]-l/2 < 1, 
This will be satisfied as long as s'^ > 0. 


( 2 . 22 ) 
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Of course, we can also decermine from q'^ and 

f^ - {l - Mq^ [4(q^)3 + (q» >3)-l (2.23) 

The problem Is that there Is no guarantee that (2.21) and (2.21) will agree. 
We can snow that the mixed triple moments s'^q' and s’q'^ provide two more 
independent formulae for fj^. With only the first and second moments ws had 
too many unknowns and not enough equations. Now with the addition of the 
third moments we have too many equations and not enough unknowns. We have 
to introduce some additonal unknowns by revising our fundamental assumptions, 
and admitting that the variances are not entirely due to a mixture of just 
two types of air. 

This we could do by Introducing.', Type 3 air, with properties ( 33 , q 3 ) and 
fraction However, this approach i’s not consistent with the observed bi- 

modal distribution of properties at cloua top. A more promising approach 
is to allow continuous if highly biraodal distributions of s and q. 
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3* Th« Coatlnuous Model 


Oc Poor qualitv 


Again we consider a oasis of ale conCalnlng two populations of parcels » 
but this. time we assume that within each population s and q have Joint noroal 
distributions about their mean values within the population (Fig. 6), l.e. 


P(*t<l) ■ f|G(s - s^yOg^t q - ^ j *"’qi*^8q|^) 


+ (I - fj) G (s - *:>*^S 2 * ^ ~ *^2****l2*'^*^2^ * 

where G is the Joint-normal distribution: 

G (x,a^,y,Oy,rjjy) = {ZitOj^Oyd - 


-1 

■» 

- 2rxyxy + 3^ 

2(1 - r,2) 

. ®x°y ®y^ 




(3.1) 


(3.2) 


and rjjy = ( Ox<3y) la the correlation between x and y within each population. 

As before, we have assigned weight f^ to population 1. and weight (1 - f|) to 
population 2* Of course. 


00 00 

I / P(8.q)dsdq - 1 . (3.3) 

••00 «*00 


since s and q are continuously distributed, we refer to this as the “con- 
tinuous model." In case the two populations are the same, the model reduces to 
the gausslan model of Deardorff and Sonmerla (1977). and Nellor (1977). 

Using the properties of the Joint-normal distribution G. we can derive the 
following relations: 
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(3.4) 

q - 

£,,, +(1 

- fi>q 

2 ’ 


(3.5) 

s’ 2 


+ (Sj- 

s) 2 j + (1 - 

fj)lOs 2 ^ + ( 82 " » 

(3.6) 

qi 2 

' ^'“< 11 " 

+ (q,- 

q) 2 ) + (1 - 

fl)[aq 22 + (q^- q) 2 ] , 

(3.7) 

s'q' 


■ 7)(q^ 

- q) + (s'q’) 

J + (1 - fi)[( 82 " 7)(q2~ 

q) + (s'q')^ 


(3.8) 

- f ^(8j- 7)[3 os^^ + (Sj- 7)2) + (1 - fiKs^- 7)[30 s 2^ + (s^- 7)2] . 

(3.9) 

q’3 « fi(9j- ^)[30q^2 + (q^_ q)2] + (j . fi)(92~ 9)[30q22 + (q2“ q) 2] , 

(3.10) 

S’2q* « f^{l<Jgj 2 + 7)2](q^ - q) + 2(8^ - 7)(s’q*)^} , 


+ (1 - fj){[Og 2 ^ + ( 8 ^- s) 2 ](q^- q) + 2 ( 8 ^ - 8 )( 8 *q') 2 } , 


(3.11) 


s'q'2 - ^^{(Sj- s){crq^2 + (q^_ q)2] + 2(q^ - q)(s»q^)^} 

+ (1 - f j){(s2 - 7)[Oq22 + (q^- q)2] + 2(q^ - q)(7^)2l • 


(3.12) 


In (3.8), (3.11), and (3.12), ( 8 *q')^ and (s'q')^ are the covarlancea associated 
with populations 1 and 2, respectively. The 9 equations (3.4 - 3.12) Involve 
the 11 unknowns sj, S 2 , qi, 92. Osj^. <^ 82 ^* ®qi^» ^ 1 * ( 8 ^)l» and (T^)^. 

Two additional equations will be required for solution of the system. The^s 
additional equations essentially determine (s'q* ) ^ and (s'q ’ ) ^» 
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If any of the triple moments are non-zero, we will necessarily find 
that 8^ fi 8^ s, or q or both. This follows from the fact that, 

for Gaussian distributions, all odd moments about the mean vanish. 

From (3,4) and (3,6), we can derive 


^ - 7+ ((s’2 -7;’^)(1 - fj)/fjl 


1/2 


■ S - t(s'^ - 0g2)f^/(l - fj)l^^^ • 


Op POOR QUALITY (3,13) 


(3,14) 


where 


Os^ = + O “ f^) Os2^ • 


(3,15) 


Similarly, from (3,5) and (3,7), we find that 


- q + M [(q'2 -7q2)(l - 


1/2 


(3.16) 


q2 “ q “ M [(q*^ - <jq2)f^/(l - fj)1^^2 , (3.17) 

where 

0^2 = fjOq^2 + (1 « ^i)Oq2^ • (3.18) 

From (3.16-3,17), we recognize the rather obvious requirements 

0s2<;”7 ^ (3.19) 

7q2<^ . (3.20) 

We now use (3,4) and (3,5) to eliminate 82 nnd q 2 from (3.6) ~ (3,12): 

7^ - Os2 + (g^« 7)2f^/(l _ f^) , (3.21) 

« Oq2 + (q^- 7)2f^/(l - f^) , (3.22) 
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7 ^ - (8^- 7)(q^- q)f^/(l - fj) + (1 - f ^)(7V)2 . (3.23) 

7^ - 8)l360g2 + (8^- 8)2 (1 - 2f j)/(l - f j)2] , (3.24) 

q^ - f j(qi“ q)[360q2 + (q^- q)2 (i _ 2f ^)/(l - f ^)2l . (3.25) 

s'2q« . f^|(q^_ q)l6og2 + (g^_ S') 2( l- 2 f ^)/( 1 -f ^) 2] + 2(8^- F)6(s'q* ) } , (3.26) 

8*q'2 « f^{(s^- s’)[6aq2 + (q^- q)2( l-2f ^)/( 1-f ^)2] + 2(q^~ q)6(s*q’ ) } . (3.27) 

Here we have used the convenient definitions 


50g2 = Og^2 - Og^^ * 

(3.28) 

«<Jq^ = Oq ^2 - <,^^2 ^ 

(3.29) 

5(s’q’) = (8*q’)^ - (8*q’)^ . 

(3.30) 


Since "og 2 and 7q2 appear only In (3.21) and (3.22), respectively, the remaining 
five equations (3.23“3.27) can be used to determine the five unknowns (sj - s), 

(qi “ q)» ^1» functions of (s'q')^ and (s'q')^. Inspection 

of (3.24-27) shows that fj cannot be zero unless all of the triple moments vanish. 
We can eliminate 60 g 2 between (3.24) and (3.26), obtaining 

7^(q^- q) - 3s^V(s^- ¥) - -2f ^(s^- 7)3(q^- q)( 1 - 2fj)/(l - f^)2 

•- 6f^ (Sj - s)2 5(7^) . (3.31) 

This can simplified using (3.23): 

• 

s'^(q^- 7) " 3s*2q»(s^- s) ■ - (s^- s) S, (3.32) 

where 

S = (s^-?){F(f^) [iTqr- (1 - + 6f^5(7V)} , (3.33) 
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F(fj) s 2(1 - 2f^)/(l - fj) . OF POOR (3.34) 

Similarly, we can use (3.23), (3.25), (3.27), and (3.33) to obtain 

q'^(s^“ ¥) - 3s'q’2(q^- q) » - <q^- q)2s/(s^- s) . (3.35) 

From (3.32) and (3.35), we obtain a cubic equation for S: 


6 s'^q' + S 3[s'q*2 • s'^ + 3(s'2q’)2j 

+s'^(s’^q'^-98*2q*8'q'2)«0. (3.36) 


We can also show that 


q3 - q 26 s'q'2 + Q3[s»2qi . q*3 + 3(gfqi2)2j 


+ q*^(s’^ q'^ - 9 s’2q» s’q'2) « 0, (3.37) 

where 


Q = (q^- q) {F(f^)[s'q'' - fj(s;q»)j- (1 - f j)(s»q») 2 l + 6fi«(s'q’)}.(3.38) 


The cubic equations (3.36) and (3.37) can be solved by standard methods (e.g., 

Abraraowltz and Segun, 1968); we choose the roots which are always real. In 

this way* S and Q can be obtained directly from the known moments. 

In the remainder of this Section, we consider the special case In which 

(s’q* ) » (FOfO “ 0. For this case, we can use (3.33) and (3.38) In (3.23) 

1 2 

to show that fj satisfies 

(16 + X)fj2 > (16 + X)f^ + 4 » 0, (3.39) 
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where 

X = SQ/(i^q"*")^ (3. AO) 

Since SQ and s^~*" always have the same sign, X cannot be negative* The solu- 
tion to (3.39) is 

fj- {l + [X/(16 + X)11''2}/2 , (3.41) 

For X =» 0, we have « 1/2. From the form of (3.41), it is clear that 
0 £ fi £ 1 will be satisfied. 

How can we choose the sign in (3.41)? For very large X, f^ approaches 
either 1 or 0, depending on the choice of sign. From (3.34), we see that F 0 
for 0 < fi < 1/2 (choice of the minus), and F £ 0 for 1/2 £ fj £ 1 (choice of 
the plus). Since, by convention, (sj - s’) 0, (3.33) implies that F must be 

positive when s’q' and S have the same sign; otherwise, F must be negative. 
Therefore, when S and s'q' have the same sign, we take the minus in (3.41); 
otherwise, we take the plus. 

If S and Q both turn out to be non-zero, completing the solution is straight- 
forward. After solving for S and Q using (3.36) and (3.37), we obtain fj from 
(3.41). Next, s-sj and q-q^ are obtained from (3.33) and (3.38), respectively. 

Of course, S 2 ~^ and q2“9 are determined from (3.4) and (3.5). Finally, Og^, Oq^, 
60 g 2 , and 60 q^ are obtained from (3.21), (3.22), (3.26), and (3.27). It is then 
easy to obtain Osj^» ®S 2 ^* ®qi^* 

On the other hand, suppose that S = 0. Then from (3.40) and (3.41) we have 
= 1/2, which implies F « 0. It follows that Q - 0. Thus, S and Q are always 
zero or non-zero together; and when they are zero, fj =« 1/2. Inspection of (3.36) 
and (3.37) shows that this happens for 

7^3 7^3 _ 9 » 0. (3.42) 
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In this case, we cannot use (3*33) and (3#38) to solve for s^ - s and 
Instead, we return to (3.23) - (3.27), which are now simplified to 


- q. 


8'~q' ' (Sj- s')(qj- q") , 

(3.43) 

(2/3) 7^ - (Sj- s) 

(3.44) 

(2/3) q'3 « (q^- q) fi0q2 , 

(3.45) 

2 s'2q' « (q^- q) 60g2 , 

(3.46) 

2 s'q' 2 - (Sj- s) 60 q^ . 

(3.47) 


From (3.45) and (3 .47), we obtain 


3 s'q’2 (q^- q) « q»3 (s^- s). 

Multiplying (3.48) through by (s^- s), and using (3.43), we find that 


s^“ s * (3 s’q*2 s'q’ / q'3)l/2 ^ 


Here we have followed our convention that sj > s. Similarly, we have 


q^- q = M (3 s'^q’ s'q' / s'^)^*^^ . 
From (3.46) and (3.47), we now find that 


60e^ * 2M (s'^q' s'^ / s'q')^^^ 


6on^’ • 2(s'q'^ q'^ / s'q')^^^ . 


If we multiply (3.49) and (3.50) together, we get 


(Sj- s)(q^- q) - M s'q' [9 s'^q' s'q'^ / (s'^ q'3)]l>'2 . 


(3.48) 


(3.49) 


(3.50) 


(3.51) 

(3.52) 


(3.53) 
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use of (3.42) now allows us to recover (3.43). 

Obviously, each quantity that appears under a square root in (3.49)-(3*52) 
must be non-negative. There is no guarantee that this will be the case; if it 
is not, then the given (or predicted) moments simply cannot be described by the 
model. This must be accepted as a possible outcome. If such a failure of the 
model occurs, it is analogous to the failure of a higher-order closure model 
chat predicts an unrealizable set of moments. 


4* Discussion and conclusion 


We have solved the continous model for the case (s'q')^ " 0» Is this 
sufficient? Such a solution will have Sjr‘s 2 andq^j‘q 2 whenever s’q’ ^ 0; 
otherwise (3.8) cannot be satisfied. Thus we will almost always have two distinct 
populations. However, the joint distribution of s and q need not necessarily 
have two maxima. Perhaps the additional generality provided by non-zero 

■ "I" 

values of (s'q* ) and (s'q * >2 is not needed. Further study is needed to resolve 
this issue. 

To apply the theory in a higher order closure model, it is necessary to 
consider the joint distributions of s and q with the vertical velocity w, 
in order to determine the buoyancy flux. Can we assume that (s,w) and (q,w) 
are jointly blmodal? This would imply that at a given level the thermodynamic 
properties of rising and sinking air have been determined by different processes— 
a very plausible idea. Further studies are in progress. 
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Figure 1: 


Figure 2: 


Figure 3: 
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Figure 5: 


Figure 6 


Schenatlc diagram illuatratlng the Irregularitlea In cloud top 
height . 

Observed Joint probability distributions of potential temperature 
and specific humidity for an aircraft flight-leg during AMTEX, as 
reported by Mahrt and Faumier (1982). Darkest shading indicates 
> 200 observations, double-hatched 100, single-hatched > 50, and 
outlined >,25. 

It 

ScheoMtic diagram illustrating the spatial distributions of Type 1 
and Type 2 air in the discrete model. 

Sketch showing the assumed joint probability distributions of 
liquid water static energy and total mixing ratio for the discrete 
model. 

Diagram used to relate changes in s*'^ to changes in fp for the 
discrete model. Along the dashed line, f^ ^ 1/2. 

Sketch ehowlng the assumed joint probability distributions of liquid 
water static energy and total mixing ratio for the continuous model. 
For Illustrative purposes, it Is assumed that the cool-wet population 
has an Internal correlation consistent with the Clauslus-Clapeyron 
relation, and relatively large standard devlationiij, while the warm- 
dry population has no internal correlatloii, and relatively small 
standard deviations. 
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